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SUMMARY 


The  effect  of  reactant  consumption  in  nonsteady  explosion  theory 
is  investigated  by  an  exact  numerical  method  and  an  approximate  analytics 
treatment.  Results  are  presented  for  the  critical  conditions  and  the  igni¬ 
tion  delay  for  supercritical  conditions  The  approximate  analysis  gives 
ignition  delays  in  good  agreement  with  the  numerical  solution,  but  the 
critical  conditions  are  not  obtained  as  accurately  as  in  an  alternative 
analysis  due  to  Thomas  A  comparison  of  temperature-time  histones 
computed  by  the  exact  numerical  method  with  experiments  reported  by 
Akita  showed  satisfactory  agreement  and  indicated  that  reactant  con¬ 


sumption  was  significant  for  wood. 
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TJl®  bafi'tc  partial  differential  equations  for  a  homogeneous  coTlh- 
btlBtldH  process  with  an  Arrhenius  type  reaction  rate  term  are 
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The  nonsteady  theory  is  an  approximate  formulation  which  replaces  the 
partial  differential  equations  by  ordinary  differential  equations  giving  the 
"average"  temperature  and  composition  as  functions  of  time.  The  space 
variation  is  eliminated  by  replacing  the  Laplacian,  which  represents  the 
heat  conduction,  by  an  equivalent  heat  loss  term  proportional  to  the 
difference  betY/een  the  "average"  temperature  and  the  ambient  temperature 
The  physical  significance  of  the  average  quantities  and  the  relation  of  the 
heat  loss  coefficient  in  the  nonsteady  theory  to  the  physical  boundary  con¬ 
dition  on  (1)  has  been  discussed  in  recent  papers  by  Gray  and  Harper^ 

* 

and  by  Thorrfes'^  » 

In  terms  of  dimensionless  variables  (defined  in  thd  List  of  Symbols) 

the  ordinary  differential  equations  are 
*  * 
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subject  to  the  conditions  0=0  and  \  =  1  at  T  =  0.  The  form  of  the 
equation  is  independent  of  geometry  which  only  eaters  in  the  determination 

of  the  lieat  transfer  coefficient,'  A 

The  object  of  this  paper  is  to  present  approximate  analytical 
and  numerical  solutions  of  (Z)  While  there  have  been  a  number  of  . 
previous  investigations,  simplifications  have  been  made,  eepecially  of 
(Z).  While  these  approximations  arc  justified  for  gaseous  reactions, 
which  were  the  motivation  for  the  earlier  work,  Thomas'^  has  pointed  out 
that  for  solids',  such  as  wood,  where 'B  is  sm'aller  by  an  order  of  magni-  . 
tude,  it  is  important  to  take  the'reactant  consumption  into  account. 

It  is  well' known  that  '(Z)  has  two  types  of  solutions  depending  on  the 
value  of  A/  6.  For  large  values  of  A/6,  the  heat  generated  by  the  reaction 
can  be  conducted  away,  and  there  is  only  a  small  temperature  rise,  the 
maximum  value  of  0  being  about  1.  This  type  of  solution  will  be  referred 
to  as  stable  or  subcritical.  For  small  values  of  A/6,  unstable  or  super¬ 
critical  solutions  are  obtained  The  heat  cannot  be  conducted  away 
rapidly  enough,  and.-,0  approaches  the  adiabatic  value  b.  The  transition  is_ 
quite  sharp  for  large  values  of  B  The  important  problems  are: 

.  3 

(1)  To  determine  the  critical  value  of  A/6  dividing  the  two 
types  of  solution,  and 


iZ)  To  determine  the  time  at  which  the  rapid  t»mp(}reij.ure 
rise  begins  in  the  supercritical  case.  Conventionally, 


s 


o  this  ignition  delay  is  defined  as  the  time  at  which 
o  Q  -  Z  For  in'aterials  such  as  wood  with  a  small  heat 

of  reaction,  there  is  no  explosion  in  the  literal  sense 
but  the  phenomenon  is  similar  in  principle  to  the  thermal 
explosion  of  gases. 


4 


0 

* 


d 


4 


^AL^fSIS 


Both  the  approjfimflte  Analytical  solution  and  the  ntimej;ical 
solution  take  advantage  of  the  absence  of  t  from  (2).  Multiplying  (2b) 
by  n\^^  ■  ^  and  dividing  by  (2a)  gives  the  differential  equation 


e 

dh^  n  \2n-lgl  +  e0 

IF  ~  ■  B  e 


subject  to  X,  =  1  at  0  =  0,  After  solving  this  equation  either  analytically 
or  numerically,  the  solution  can  be  expressed  as  a  function  of  t  by  the 
auxiliary  relation 
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X"e 


1  +e0 
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The  use  of  0  as  the  independent  variable  involves  some  problems  as  both 
X  and  T  are  double  valued  functions  of  0; --howeYOr ,  these  are  not  significant^ 
in  the  specific  problems  under  consideration,  and  ibwas  found  easier  to 
work  with  (3)  than  with  the  alternate  form 


d0 

dX 


0 

A  0  '1  +  00 

T - e 

■6  ),n 


(5)  . 


Ncmerfcal  SoliltlQn 


C< 


In  this  section, only  a  first  order  equation  (n  =  1)  is  ffOnsiderecJ, 

* 

and  e  is  neglected,  With  these  simplTfications,  (3)  can  be  written  as  an 
integral  equation 


1  -  B 
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(6) 


This  integral  form  is  only  valid  while  9  increases  to  its  maximum  value 
However,  no  difficulty  results  from  this  limitation  In  subcritical  cases 
an  imaginary  value  of  X  is  obtained  when  the  calculation  is  carried  past 
the  maximum  value  of  0  In  the  supercritical  cases,  the  calculation  is. 
only  carried  out  to  0  =  2  which  is  well  below  the  maximum  value. 

The  numerical  solution  is  based  on  the  relations 


2X(0)  = 


f-0e-0  +  F,0 


4^0e-^  +  F(0  -  A0) 


.  4  A  0e'®  F(0  -  A0) 
2B/  6 


1/2 


(7a) 


where 


A0 

F(0  -  A0)  =  X(0  -  A0)  -  j- 


_ M0  -  A0) _ 

X(0  -  AO)  -j  (0  -  A0)e-(®  - 


(7  b) 


This  pfocedure  is  equ5x»alertt  to  the  veil  ■  lirtown  tnodifled  Eulej* 
rfiethod^  of  solving  an  ordinary  differential  equation  The  iterative  pro¬ 
cess  was  eliminated  by  the  assumption  n  =  1  which  made  the  nonlinear 
algebraic  relation  a  quadratic  which  was  solved  to  give  an  e^fplicit 
formula  for  \{d)  in  terms  of  ^(0  -  A0)  and  0.  The  method  is  also  appli¬ 
cable  with  c  taken  into  account  The  simplifying  a'ssumption  was  made 
because  e  is  in  fact  small  and  the  approximate  analytical  solutions,  with 
which  the  numerical  solution  will  be  compared,  neglect  e 

It  can  be  seen  that  the  numerical  'Solution  of  the  basic  equations 
only  solves  the  first  problem  dfetermining  the  critical  value  of  A/6  for  a 
given  value  of  B,  indirectly  The  critical  value  of  A/6  for  a  given  value 
of  B  can  be  bounded  as  closely  as  desir.ed  by  means  of  numerical  solutions,, 
but  no  explicit  formula  for  the  critical  curve  can  be  obtained  It  is, 
therefore,  desirable  to  consider  approx-imate  analytical  solutions. 
Approximate  Analytical  Solution 

The  approximate  analytical  solution  of  (6)  is  based  on  the  method 
of  successive  approximations.  If. an  approximate  expression  for  X.(0) 
is 'put  into  the  integral  on  the  right-hand  side  of  (6),  an  improved  approxi¬ 
mation  is  obtained  by  evaluating  the  integral  The  first  approximation  will 


where  q  iS  a  parameter  to  “be  determined  by  matching  the  first  and  second 
approximations  at  a  reasonable  value  of  0.  Using  the  initial  slope  gives 
q  =  1 ,  but  the  higher  value  obtained  b/  a  matching  procedure  appears 

desirable. 

It  should  also  be  noted  that  within  the  accuracy  of  the  approximation 


(9) 


Therefore,  the  approximate  analysis  can  be  extended  to  an  nth  order 
reaction  by  replacing  B  by  B/n  except  when  B  is  small. 

Examination  shows  that  with  K  given  by  (8),  or  even  with  X  =  1, 
the  integral,  in  (6)  cannot  be  evaluated  in  closed  form  For  X  =•  1  (no 
reactant  consumption)  Gray  and  Harper^  used  the  approximation 

=  T+  .720  +  0^  •  (10) 

This  does  no  allow  integration  in  closed  form  when  (8)  is  used,  therefore, 
an  analogous  approximation 


1  +  .7180  ■ 

1  -  3680 

0 


*  Q 

•■■The  constants  7180- e- Zand  368  ^  1 /e  were  selected  to  give  the 
correct  result  for  B  — ►  oo 


(11)-^= 


o 


n 

■was  introduced  irt  otder  to  ctjtaic  ait  tmegrftbT<?  form.  F?g.  1,  winch 
compares  the  two  approximations  within  the  range  0^  6  2  shows  that 

the/  are  of  comparable  accuracy. 

With  these  approximations,  the  relations  - 
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d0 
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are  obtaine.d  These  integrals  can  be  evaluaj;ed  in  closed  form  by  well- 
known  methods)  though  the,  resulting  expressions  are  rather  long. 

The  Critical  Case 


The  critical  case  corresponds  to  the  denominator  in  (12),  being 
a  perfect  square.  Algebraic  manipulations  gives  the  relation 


f)  =(l. 454  -  1)  +y(l. 454  -1)8. ( 

cr  ^  , 


13) 


a 


i 


TO 


for  tKe  critical  vaTwe  of  AfS  However,  the  expression  does  not  give  S 

definite  result  until  q  is  evaltiated.  Introducing  (13)  for  A/6  into  (12a) 
gives 


X 


718  -  I)  (PO) 
(1  -  pe)^ 


718P-^  ^  (P0)^ 


d(P0) 


(14) 


where 


Physically,  P‘^  represents  th-e  maximum  value  of  0  in  the  critical  case. 
Therefore,  using  a  constant  value  of  P0  as  the  matching  point  appears 
plausible.  Fig.  2  shows  the  variation  of  q  with  B  for  several  choices 
of  values  of  P0  at  which  (8)  and  (14)  match.  It  can  be  seen  that  q  is 
constant  until  B  becomes  relatively  small,  A  discussion  of  the  resulting 
critical  curve  will  be  given  below. 

Supercritical  Cas£ 

In  the  supercritical  case,  the  value  of  A/'6  is  arbitrary,  provided 

L  s> 

it  is  below  the  critical  value.  The  problem  is  to  compute-' the  value  of 
6t  for  0  =  2  by  evaluation  of  (12b).  It  is  reasonable  to  evaluate  q  by 


IZ 


matcTiing  at  &  -  1,  the  midpoint.  The  resulting  \raltieB  of  q  for  eeveraT  • 
values  of  A/6  are  shown  in  Fig  3  which  plots  q  versus  B/q  It  should 

be  noted  that  the  curves  end  abruptly  Corresponding  to  (13),  the  relation 


can  be  obtained,  and  for  a  given  A/6  there  are  no  supercritical  points 
for  values  of  B/q  smaller  than  those  implied  by  (15). 

.  I 


14 


SPECIFIC  RESULTS  •  ■ 

The  integral  equation  (6)  was  solved  numerically  for  a  number  of 
sets  of  A/6  and  B  using  AO  =  0.1.  Some  check  calculations  with  smaller 
values  of  A  0  (0.  05  and  0. 01)  indicated  that  the  step  size  was  sufficiently 
small.  The  results  of  these  numerical  calculations  will  be  compared  with 
the  approximate  analytical  calculations  of  the  critical  conditions  and 
ignition  delay  obtained  in  the  previous  section  and  those  reported  by 
Thomas-^.  Also,  a  comparison  with  experimental  temperature  histories 
obtained  .by  Akita^  will  be  presented. 


Critical  Conditions 


6cr(B) 


•Fig.  4  compares  the  value  of  ^ 


cr 


00 


=  e  -T  -T-  obtained 


cr 


from  (13),  the-  value  of  q  being  .obtained  by  matching  at  P0  =  0.5,  with 
Thomas' -expression^ 
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and  with  Frank-Kamenetskii' s  expression. 
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FIGURE  4.  EFFECT  OF  REACTANT  CONSUMATiON  ON 
CRITICAL  CONDITIONS  FOR  A  THERMAL  EXPLOSION 


The  limits  of  the  stability  ctli*ve  foUttd  b7  the  nitiTleflCal  ealCuTatlOilS  for 


particular  values  of  B  are  indicated  by  |  .  While  the  limits  were  not 
determined  closely  except  for  B  =  176,  it  is  clear  that  Thomas'  result 
is  definitely  superior  to  the  others  except  for  large  B  when  the  effect  of 


reactant  consumption  is  small. 


Ignition  Delay 


The  ignition  delay  T.^  is  defined  somewhat  arbitrarily  as  the  time 

required  for  0  to  reach  Z.  Fig.  B  shows  as  a  function  of  A/6  for 

B  =  Z9. 4  and  B  =  176  as  obtained  by 

(1)  Numerical  solution  of  (6)  combined  with  numerical 
evaluation  of  (4) 

(Z)  The  approximate  analysis  based  on  (IZb)  with  q 
determined  by  matching  at  0  =  1 

(3)  Thomas’  approximate  analysis,  and 

'  (4)  Numerical  integrations  carried  out  by  Rice,  Allen  and 
7 

Campbell  in  which,  however,  reactant  consumption  is 
not.  taken  into  account  rigorously. 


It  can  be  se’en  that  in  this  case  our  approximate  analysis  is  in 


good  agreement  with  our  numerical  calculations  .The  approximate 
analysis  of  Thomas  and  the  numerical  calculations  of  Rice,  Allen  and 
Campbell  are  somewhat  further  off,  but  in  the  opposite  direction,  than 


our  approximate  treatment. 

Comparison  with  Experimental  Temperature  Histories 

4 

Kinbara  and  Akita  have  reported  two  temperature  histories  for 

O 

sawdust  spheres  (contained  in  w'lre  gauze)  and  compared  the  data  with 


eoldiloTis  of  a  linc?ai*t nation  of  (la)  Th&  ag<e(»mf.'Tit  ^^^as  fiualitativel^ 
satisfactory,  but  the  ignition  delay  was  underegtimated  somewhat.  It 
was  suggested  that  this  wag  due  to  the  neglect  of  reactant  consumption. 

In  Fig  .  6a  and  6b,  the  experimental  results  are  shown  along  with  Kinbara 
and  Akita's  calculated  curves  and  numerical  solutions  for  B  =  10  and  two 
lower  values  of  B.  It  can  be  seen  that  a  value  of  B  of  about  5  leads  to  a 
marked  improvement  in  the  agreement  between  theory  and  experiments. 
Unfortunately,  Kinbara  and  Akita  do  not  give  the  value  of  Q  (the  heat  of 
decomposition  of  their  sawdust),  but  B  =  5  does  correspond  to  a- reasonable 


value . 


TEMPERATURE  RISE 
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EXPERIMENTA 


2  4  6  8  10  12  14  16  18  20 

,  TIWE(MINUTES) 

6a-3.7cm  SPHERE  AT  I96°C;  A/S  =.46 


FIGURE  6.  COMPARISON  OF  EXPERIMENTAL  AND 
THEORETICAL  TEMPERATURE  HISTORIES 


DlSCtJSSIOM 


Tlie  mathematical  analysis  of  the  differential  equations  of  non- 
sLeady  e^cplosion  theory  has  been  greatly  simplified  by  using  the  tempera¬ 
ture  instead  of  the  time  as  the  independent  variable.  In  particular,  for 
a  first  order  reaction,  the  numerical  solution  is  greatly  simplified  because 
■  a  nonlinear  relation  reduces  to  a  quadratic  which  can  be  solved  explicitly.’ 

In  the  approximate  analytical  treatment  of  the  equations,  two 
approximations  are  involved.  It  is  believed  that  the  principal  source 
of  error  is  the  use  of  the  linear  relation,  (8),  between  composition  and 
temperature  as  a  first  approximation  to  the  solution  rather  than  the  use  ■ 
of  the  rational  fraction,  (11),  to  approximate  the  exponential.  In  particula 
it  could  be  anticipated  that  the  linear  approximation  between  composition 
and  temperature  would  not  be  accurate  near  a  temperature  maximum. 

This  probably  explains  the  error  in  the  determination  of  the  critical  con¬ 
ditions  where  a  temperature  maximum  is  involved  and  the  relative 
accuracy  of  the  determination  of  ignition  delay  for  supercritical  conditions 

a 

where  no  maximum  is  involved 

It  is,  of  course,  true  that,  since  there  is  an  adjustable  parameter, 
q,  involved  in  the  calculation,  the  determination  of  the  critical  condition 
could  probably  be  improved  by  some  other  method  of  evaluating  q.  The 

O 

method  used  in  this  paper  of  matching  at  a  fixed  yalue  of  P(3  gave  a  value 


21 


of  q  whicTi  waff  essentially  independent  of  B  An  examinatlcn  of  Thomas* 

expression  indicated  that. to  improve  our  result,  q  should  increase  as 
1/3 

B  ,  so  that  no  basic  improvement  can  be  expected  from  changing  the 
value  of  P0 .  While  other  methods  of  evaluating  q  are  possible  (such  as 
matching  at  a  fixed  value  of'0  rather  than  PO),  it  hardly  seems  worthwhile 
to  attempt  to  patch  up  the  analysis  in  this  way,  particularly  since  the 
exact  numerical  solution  is  so  convenient. 

The  comparison  between  Akita' s  experiments  and  the  calculations 
indicate  that  the  agreement  is  improved  by  taking  reactant  consumption 
into  account,  but, as  Akita  did  not  specify  the  value  of  Q  for  the  sawdust, 
a  question  remains  as  to  whether  the  value  of  B  which  gives  satisfactory 
agreement  is  the  actual  value  for  the  material  though  it  appears  reasonable, 

5  f 

In  view  of  the  uncertainty  'in  the  measurements  and  in  the  value  of  the  other 

properties  involved,  the  agreement  appears  as  good  as  can  be^expected. 

•  * 

It  v/ould,  however,  be  of  some  interest  to  treat  the  problem  by  numd'rical 

u 

solution  of  the  full  pai-lial  differential  equations  (la)  and  (lb)  in  order 
to  estimate  the  effect  of  the  approximations  involved  in  nonsteady 
explosion  theory,  -A  numerical  procedure  for  solving  the  partial  differential 

Q 

equation  has  been  developed®  but  only  a  few  preliminary  calculations  have 
been  made  as  yet 

The  numerical  solutions  for  small  values  of  B«  10)  show  that 
the  transition  between  the  super-  and  subr.ritical  cases  becomes  less 


distinct  as  iB  ‘becemiBS  Smaii.  Exanii^iation  of  ft 4)  sTio'X's  tTiat  tho  niaxi- 
mum  temperature  in  the  stibcritical  case  mcreases  as  B  decreases,  while 
the  maximum  temperature  in  eupercritical  qasee  must  be  less  than  B,  so 

that  the  distinction  obviously  becomes  less  sharp  for  small  values  of  B 
In  this  region,  it  is  necessary  to  use  the  numerical  procedure  rather  than 
the  analytical  approximations  which  are  essentially  ejcpansions  in  B"^. 
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